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Abstract 

The coherence time, and thus sensitivity, of trapped atom interferometers that use non- 
degenerate gases are limited by the collisions between the atoms. An analytic model that describes 
the effects of collisions between atoms in an interferometer is developed. It is then applied to an 
interferometer using a harmonically trapped non-degenerate atomic gas that is manipulated with a 
single set of standing wave laser pulses. The model is used to find the optimal operating conditions 
of the interferometer and direct Monte-Carlo simulation of the interferometer is used to verify the 
analytic model. 
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I. INTRODUCTION 



To date cold atom interferometers have demonstrated rotation sensitivities comparable 
to ring laser and mechanical gyroscopes Several atom interferometer schemes have been 
realized, and thus far free space fountain and beam configurations, that utilize light pulses to 
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While the lack 



manipulate the atoms, have demonstrated the greatest sensitivities 
of external potential reduces systematic errors, atom interferometry in free space is limited 
by acceleration due to gravity. In particular, the precision of an interferometer is directly 
proportional to the interrogation time and in free space this time is limited by the size of 
the vacuum chamber. In the most sensitive free space atom interferometers, the atomic 
clouds travel up to 10 meters The large scale of free space interferometers limits their 
applications. 

There is currently a great effort being made to reduce the size of atom interferometers 
while simultaneously increasing their sensitivity. One straightforward way to achieve this 
goal is to develop interferometers that trap the atoms in an external potential for the duration 
of the interferometer cycle. The external potential prevents the atoms from falling due to 
gravity, and keeps the atomic gas from expanding in the vacuum chamber. As a result, the 
interferometer cycle time is not as limited by the size of the chamber. 

Several groups have built trapped atom interferometers using atomic gases that are both 
above and below the recoil temperature . To date, all interferomters that 

use gases below the recoil temperature have utilized atoms in a nearly pure Bose-Einstein 
condensate (BEG). 

If the gas is cooled below the recoil temperature and is split using a laser pulse, a large 
relative separation between each arm of the interferometer can be achieved [^. By exposing 
the atoms in each arm to a different environment, precision measurements of localized phe- 
nomena can be performed. For example, the AC stark shift in ^''Rb was recently measured 
by exposing the atoms in one arm of an interferometer to laser light |lO |. 

In some applications, such as the sensing of rotations and accelerations, cloud separation 
is not necessary and interferomters that use non-degenerate source are sufficient. These 
sources can be produced by laser cooling alone. Additionally, non-degenerate atomic gases 
have a much lower density compared to a BEC and therefore experience a weaker mean- 
field potential. The mean-field potential directly couples number uncertainty into phase 
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uncertainty via number dependent phase diffusion [8|, [111, Il2|, Il3| . One advantage of working 
with laser cooled gases is that dephasing due to number fluctuations is ameliorated. 

Besides the elimination of mean-field effects, atom interferometers that use laser cooled 
atomic gases are less sensitive to heating due to imperfections in the confining potential than 
a BEG. Atoms in a pure BEG experience no momentum changing collisions with other atoms 
in the same mode. Therefore, an interferometer that uses a pure BEG will experience little 
decoherence due to collisions. However, if a BEG is heated, atoms will leave the condensate 
and will experience an increase in the collision frequency. As a result, the decoherence rate 
will increase with temperature. On the other hand, if the interferometer uses a laser cooled 
gas with a temperature much greater than the BEG transition temperature, the density will 
decrease as the temperature increases and the decoherence rate will decrease if the gas is 
inadvertently heated. 

Several different methods for building atom interferometers using laser cooled gases have 
been developed and time-domain atom interferometers that use a single internal quan- 
tum state 15| lend themselves naturally for use with trapped atomic gases. This type of 
interferometer uses a series of optical standing waves to manipulate the external states of 
the atoms in the cloud. The interferometric cycle begins by loading an atomic cloud in a 
magneto-optical trap. The trap is switched off and the atomic cloud begins to fall due to 
gravity. At the time t = 0, the cloud is illuminated with a short pulse from the standing 
wave laser field. Shortly after the pulse, the cloud has a density modulation with a period 
of A/2, where A is the wavelength of the laser field. The density modulation then disap- 
pears because of the thermal motion of the atoms in the cloud. At the time t = T the gas 
is illuminated with a second pulse. Due to the Talbot-Lau effect, there is an echo of the 
density modulation at the times t = nT, for integers n > 2. If the atomic cloud experiences 
a non-uniform potential during the interferometer cycle, the density echos will be shifted 
relative to the initial modulation. The shift in the phase of the modulation can be deter- 
mined by reflecting a probe pulse from a single laser beam off of the echo. The phase of the 
reflected probe pulse is directly proportional to the phase shift in the density modulation. 
By interfering the reflected probe pulse with a reference beam, the interferometer's signal 
can be read. 

Recently, a trapped time domain atom interferometer was built by the group at Harvard 
Q]. This interferometer used an atom wave guide to confine the atoms in the perpendicular 
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directions while allowing them to freely propagate along the parallel direction. A series of 
standing wave laser pulses were applied to the atoms, such that the wave vectors of the 
lasers pointed along the free direction of the guide. The Harvard group demonstrated that 
it is possible to electronically move the wave guide back and fourth perpendicular to the free 
direction of the wave guide so that the arms of the interferometer enclose an area, making 
the interferometer sensitive to rotations j^. 

A major difficulty with all trapped atom interferometers that use optical pulses is that the 
residual potential along the guide causes decoherence [itI . [isl . [19I . The groups that have 
built BEC based interferometers have mitigated the decoherence by either using a double 
reflection geometry or using the classical turning points of the residual potential to reflect 
the atoms. The Harvard group has reduced the effects of the residual potential by using 
an interferometric cycles with several laser pulses 2^. Although this multi-pulse scheme 
greatly increases the coherence time of the interferometer, it also reduces the number of 
atoms participating as well as reducing the area inclosed by the interferometer. 

We are currently developing a trapped atom gyroscope that uses a laser cooled atomic 
gas and avoids decoherence due to the residual potential by using classical turning points to 
reflect the atoms. Rather than utilizing the Talbot-Lau effect, the density modulation will 
echo twice every oscillation of the atoms in the parallel direction. Figure [1] is a schematic of 
the interferometer cycle, (a) Initially, a laser cooled atomic gas is loaded into a cigar shaped 
trap. The trap in the perpendicular direction is created with the upper most horizontal wire 
plus a uniform bias field. The relatively weak trap in the parallel direction is created using 
vertical wires that are not shown, (b) At the beginning of the interferometer cycle t = 0, 
the atomic cloud is illuminated with a standing wave laser field. The atoms are accelerated 
towards the nodes of the laser field, (c) Immediately after the laser pulse, the atoms move 
towards the location of the nodes and density modulation appears across the cloud, (d) The 
density modulation disappears due to the thermal motion of the atoms. Simultaneously 
the trap is moved downwards by cycling the current in the wires. If the interferometer is 
rotating about the plane of the paper with frequency Q, the Coriolis force will accelerate 
the cloud in the parallel direction, (e) The cycling of the currents in the wires is timed 
so that the trap is above the bottom wire at half a period of the parallel trap t = T/2, 
where T is the trap period in the parallel direction. Near t = T/2 there in an echo of the 
density modulation across the cloud, (f) The trap is moved upward. The Coriolis force 
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decelerates the cloud resulting in a displacement of the cloud that is directly proportional to 
the rotation frequency of the interferometer, (g) At the time t = T, the trap returns to the 
top wire and there is a second echo of the density modulation. The modulation is shifted 
due to the rotation of the interferometer. The cycle (c) through (g) can be repeated many 
times. Since the oscillating Coriolis force is resonant with the parallel trap frequency, the 
shift in the displacement of the cloud will increase after each cycle, (h) After n cycles, the 
displacement of the cloud is precisely measured by reflecting a probe beam off of the density 
modulation and interfering the reflected light with a reference beam. 

The probe pulse only interacts strongly with the cloud when there is a density modulation 
across the cloud. As a result, the probe pulse can be longer than the duration of the 
modulation echo. Small fluctuations in the trap frequency can be measured simultaneously 
with the phase shift and using thermal atoms avoids the critical timing needed when a BEC 
is used fl fl Q. It „,a. p..ove possible to .easu.e the i.terfete„ee signal .no.e than 
once in any given experiment. As a result, it might be possible to split the atoms once, and 
measure the rotation frequency several times as the cloud oscillates in the trap. 

During interferometer cycle, collisions between the trapped atoms will bring the gas back 
to equilibrium, causing a reduction in the amplitude of the density modulation. Thus, the 
amplitude of the reflected probe pulse will degrade with time. The upper limit on the 
interferometer cycle time and the devices sensitivity can be determined by analyzing the 
effects collisions between the atoms in the trap. 

In this paper, we present a theoretical model for our interferometer. In Sec. [Tll we present 
an analytic model for the amplitude of the reflected probe pulse including the rotation of 
the interferometer and the effects of collisions between the atoms in the gas. In Sec. IIIII 
the analytic model will be used to estimate the minimum value of the trap frequency in 
the perpendicular direction, the optimal temperature of the gas, and the optimal number of 
atoms to use in our upcoming experiment and we compare the results of our analytic model 
with a Direct simulation Monte-Carlo code. Finally Sec. IIVI conclusions will be presented. 
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II. FORMULATION OF THE PROBLEM 



The dynamics of a dilute atomic gas a 



by the quantum Boltzmann equation 



)ove BEC phase transition temperature is governed 
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|.4p.„,.1 = /(p). (1) 

where p is the single particle density operator, -ffe// is the effective single particle Hamilto- 
nian, and /(p) is the collision integral. The effective Hamiltonian for an atom in a rotating 
trap and standing wave laser field is 

2 

Heff = ^ + V{r) + hxcos{2ki ■ r) - H ■ (r X p) + 2f/on(r), (2) 

where m is the atomic mass, V is the external trapping potential, x characterizes the strength 
of the standing wave laser, ki is the wave vector of the laser field, and f2 is the vector that 
points along the axes of rotation with the magnitude of the angular rotation frequency. The 
final term in Eq. ([2]) is the mean-field potential where Uq = Anh'^as/'m characterizes strength 
of the atom-atom interactions, is the s-wave scattering length, and n is the number density 
of the atomic gas. Note that the mean field potential for a non-condensed gas is a factor of 
two larger than for a BEC with the same density. 

It is convenient to recast the single particle density operator in the Wigner function 
representation which is defined as 

/(r,p) = ^ / d'r'{r-r'\p\r+r')e'^^'-^/\ (3) 

where \r) are the eigenstates of the coordinate operator. The Wigner function can be 
interpreted as the probability density of finding an atom at the coordinate r with momentum 
P- 

It will be assumed that the standing wave laser pulse is in the Kapitza-Dirac regime. It 
is sufficiently short that both the free evolution of the gas and the collision integral may be 
neglected, i.e. the atoms do not move and experience no collisions while the laser beams are 
on. The pulse is in this regime when 

rp«-,rp«-, (4) 

where Tp is the length of the pulse, A is the wavelength of the laser beams, v is the average 
speed of the atoms in the gas, and z/ is the average collision frequency. When Eq. (jl]) is 



fulfilled the dynamics of the atomic gas can be separated into two parts: the dynamics when 
the laser beams are on and the dynamics when the laser beams are off. 

In what follows the dimensionless coordinate r' = 2kiX, the dimensionless momentum 
p' = p/2hki and the dimensionless time t' = t/to where to = fn/Ahkf will be used. For 
'^^Rb, the characteristic time is to = 5.3 /is. Substituting Eq. and into Eq. ([T]) the 
dimensionless equation of motion for the Wigner function /, when the laser beams are on, 
is 

= Xsin(2x)[/(r,p-A:,/2)-/(r,p + V2)], (5) 

where fcj = ki/ki is the direction of the standing wave laser field, x' = toX is the dimensionless 
laser strength, and all of the primes have been dropped. Similarly, the dimensionless equation 
of motion for the Wigner function /, when the laser beams are off, is 

d ^ df dVeff Of 

f(r,p,t) = —p- — -\ IL . — 

5t ' ' dr Or dp 

- n.(rx^+px^) (6) 

where Vt' = to^, y'o = 2kiyo and once again all the primes have been dropped. Our short 
term goal is to measure the rotation of the Earth. For ^''Rb the rotation frequency of the 
Earth in our dimensionless units is fi^; = 4 x 10~^°. When the gas is in thermodynamic 
equilibrium, the dimensionless temperature is T' = T/ATpt, where Tr = h'^kf/mkB is the 
one photon recoil temperature. For ^^Rb, the recoil temperature is = 350 nK. 
The dimensionless effective potential is 

Veff = V + 2gn, (7) 

where V = toV/h is the dimensionless trapping potential, n' = 8kfn is the dimensionless 
density, and g = Snagki is the dimensionless mean-field strength. For ^''Rb the dimensionless 
mean-field strength is ~ 1. 

The length scale L of density changes in a magnetically trapped atomic gas is typically 
much larger than the atoms s-wave scattering length a^, i.e. a^/L ^ 1. In this limit, the 
collision integral becomes independent of the potential. Since the atomic gas is above the 
BEG transition temperature, no single quantum state has a macroscopic population and 
Bose enhanced scattering can be neglected. The dimensionless collision integral can be 
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approximated with the classical collision integral [22 




d^P3dn\p3-p\ [/(r,pi)/(r,p2) - f{r,p3)f{r,p)] , 



(8) 



where a = 32'n'alkf is the dimensionless collision cross section. For '^^Rb the dimensionless 
scattering cross section is a = 0.2. 

Before the laser pulse is applied, the gas is in thermodynamic equilibrium and it will be 
assumed that the laser pulse is sufficiently weak that the gas is always close to equilibrium. 
The Wigner function / can be written as 



where /o is the equilibrium Wigner function and 6f is the disturbance caused by the laser 
pulse. When the disturbance is much smaller than the equilibrium \6f\ <^ |/|, Eq. (IHl) can 
be approximated as [3] 




Icon = -r / d^Psdnip, -p\ [2/o(r,pi)<5/(r,p2) - fo{r,Ps)6f{r,p) - fo{r,p)6f{r,ps)] . (10) 



Equation ffTOj) is the sum of three terms, each with a simple physical interpretation. The 
first term 2fo{r,pi)6f{r,p2) is proportional to the rate that an atom in /o scatters with an 
atom in 6f and one of the atoms scatters into the momentum state p. The second term 
fo(r,p3)6f(r,p) is proportional to the rate that atoms scatter out of 6f because of collisions 
with atoms in /q. The final term fo{f,p)Sf{f,P3) is the inverse of the second process. 

Only atoms in the disturbance contribute to the interference signal. Therefore, once an 
atom scatters out of the disturbance it no longer contributes to the interference signal. As a 
result only the second term in Eq ( ITOl) contributes to the loss of the interference signal and 
the collision integral Eq. fllOl) becomes 



/ = /o + Sf, 



(9) 





hoii = i'{r,p){fo{r,p) - f{r,p)), 



(11) 



where the collision frequency u is given by the integral 




(12) 



Substituting the equilibrium distribution 



/o 



1 




(13) 



(27rT)3/2 



n 
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into Eq. (fT2|) . the collision frequency Eq. (fT2|) becomes 



v=-p{r)aT'/'K{\p\/V2T), (14) 

TT 



where 



ir(0 = j dr]d9ri'^ sm9^/C'^ + 7]^- 2^r] cos ^e"'?'. (15) 

The integral -f^(0 can be explicitly evaluated in terms of error functions. To remove the 
dependnce of u on the coordinate r and momentum p, Eq. (fT5|) will be replaced by its value 
a zero argument K = 1 and the density n will be replaced by the averaged density of a 
gas thermodynamic equilibrium in a harmonic potential. The collision frequency Eq. (I12p 
becomes 

^^=(2^0^^' 

where u = (cuncj^)^/^ is the geometric average of the trap frequencies. In Sec. IIIII it will 
be demonstrated that using Eq. flTBl) for the collision frequency yields accurate results when 
compared to a more complete description of the atomic collisions. 

For the rest of this paper, we will limit the discussion to the case of a cigar shaped 
harmonic potential. The dimensionless trapping potential becomes 

V=l{u;lx' + col[iy-yoY + z']}, (17) 

where uj'j_ = Iquo^ and uj'^^ = tQUj\\. If the trap has frequencies uj\\ = 27r x 3 Hz and ujj_ = 
2-71 X 300 Hz, the dimensionless trap frequencies for ^'^Rb are uj'^^ = 10~^ and uj'j_ = 10~^. 
When the gas is close to thermodynamic equilibrium in a harmonic trap the density is 

where is the number of atoms in the trap and V is the potential. Using Eq. ([7]) and f[T5]) 
the effective potential can be expanded to fourth order as 

V^ff -[^- J^^2 ) y + (2vr)3/2TV2^ • (19) 
The lowest order mean field contribution to the potential causes a small reduction of the 
trap frequency. Therefore, the oscillation period of atoms in the trap is weakly dependent 
on the number of trapped atoms. The next higher order contribution is a weak quartic 
contribution to the potential. This, and all higher order terms, can be neglected when 

2guj^N 



(27r)3/2T7/2 



< 1. (20) 



For example if the trap contains 7x 10^ ^^Rb atoms, in a trap with frequencies uj\\ = 27rx 3 Hz 
and u± = 27T X 300 Hz, and a temperature of 40 fiK (T = 30), the left hand side of Eq. fl20|) 
is about 6 x 10"^. In this case, the quantric contribution to the potential can be neglected. 

The analysis of the operation of the interferometer will be limited to the case where the 
splitting and read lasers beams are aligned with the weak axis of the harmonic potential, 
which will be chosen to be the x direction. For definiteness, the rotation of the interferometer 
will be in the z direction and the trap will be moved in the y direction. An atomic cloud at 
temperature T remains in equilibrium if the center of the trap is translated adiabatically. 
The trajectory of the moving trap yo{t) is adiabatic when 

^ « v^-. (21) 

where u± is the trap frequency in the y direction. 

When Eq. ( |2T1) is fulled, the equations of motion Eq. (JSj) and ([6]) can be recast in a 
one-dimensional form. The Wigner function is written as the product f{x,r±,p,p±) = 
f{x,p)F{r±,pj_), where F{rj_,p±) is the equilibrium distribution in the perpendicular di- 
rection, normalized to one, and f{x,p) is the non-equilibrium distribution in the parallel 
direction, normalized to the number of atoms in the trap. When the laser beams are on, the 
one-dimensional equation of motion for the Wigner function is 

df 

= xsin(x) - 1/2) - fix,p+ 1/2)] , (22) 

The solution of Eq. (!22l) can be written in terms of Bessel functions of the first kind J^, 

/(x,p,t) = 5^(-z)'j,(H)J,+,(H)e^('+2^)-/o(x,p-//2), (23) 

Ik 

where the sum over k and / runs from — oo to oo, 5 = J dt'xit') is the strength of the laser 
pulse, and 

(24) 



2T 

is the equilibrium Wigner function at temperature T. In general Eq. fl23|) can be negative, 
because the resulting gas is in a non-classical state. However, for high temperatures T 3> 
1 the negative parts of the Wigner function are negligible, and the gas may be treated 
classically. 

After the laser pulse is applied, the optical field is turned off and the one-dimensional 
equation of motion for the Wigner function becomes 
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/ = K/o-/), (25) 



where the colhsion frequency u is given by Eq. (fT6!) . The left hand side of Eq. (!25l) can be 
greatly simplified by introducing the new coordinates 

P /"* 

x' = X cos LUut sintuiit + i7 / dTyoir) cos luiit 

^11 J 
p' = C(J||Xsinc(j||t + pcosc(J||t + i7ci;|| J dryQ^r) sin cu^^t 
t' = t. (26) 

In this new coordinate system, Eq. 0251) becomes 

^/ = K/o-/), (27) 

which has the general solution 

fit) = f{0)e"'' + fo, (28) 

where /(O) is the initial Wigner function given by Eq. (123|1 . /o is the equilibrium Wigner 
function given by Eq. flM|) . and u is given by Eq. f|T6|) . 

To read out the accumulated phase, the atomic cloud is illuminated with a single off 
resonate laser beam. The light that is back scattered off of the cloud is mixed with a 



reference beam 
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171]. By measuring the interference intensity, the amplitude of the 



scattered light can be determined. Using the Born approximation, it can be shown that the 
amplitude of the back scattered light is proportional to [l3] 

S = j dxdpe'''f{x,p), (29) 

where / is the one-dimensional Wigner function. The quantity S will be referred to as the 
interference signal of the interferometer. 

In the new coordinate system Eq. ( !26|) . the signal becomes 



S = e-'^ j c/a;c/pe*(^''=°'"*+P/"''''^*V (30) 

where 

(p = Q^osu!t J dt'yo(t') cos tot' + sin tot J dt'yo{t') sintot'l^ (31) 

11 



is the phase shift due to the rotation of the interferometer. 
Substituting Eq. ([28D into Eq. ([30]) yields 

I nfe J 

The interference signal S is only nonzero when t ~ 2mT/uj, where n is an integer or half 
integer. Expanding Eq. fl32|) near these points and taking the limit where T ^ uj"^ yields 

00 

S = -2iVe-^^-^"'/2-^^t J2 JkJk+i sm{{k + l)r/2) (33) 

k=0 

where t = t — 2n'K /uj\\. Equation (!33|) along with Eq. (fT6|) and (13T!) are the main analytical 
results of this paper. 



III. DISCUSSION 

Equation (l33l) will now be analyzed and optimal operating conditions for trapped thermal 
atom interferometers will be found. Additionally to confirm the results of the analytic model 
we will use a direct simulation Monte-Carlo model (DSMC) of the interferometer 24 1. 

Direct simulation Monte-Carlo is accurate because the equation of motion after the laser 
pulse Eq. ([6]) is equivalent to the classical Boltzmann equation. Since the effect of the 
standing wave laser pulse on the cloud is non-classical (Eq. (122!) ), DSMC can only be used 
to model the dynamics when the standing wave laser beams are off. To account for the laser 
pulse, the initial conditions for the DSMC model was set by Eq. (1231) . In regions where the 
initial Wigner function is negative / < 0, the classical distribution of atoms was set to zero. 
This is valid when the temperature is much larger than the two-photon recoil temperature 

r> 1. 

For definitiveness, we specialize to the case where the trap is moved back and fourth 
according to 



Z/o = ^cos(u;||t), (34) 

where d is the dimensionless distance that the atomic cloud is displaced in the ?/-direction. 
Our chip will displace the trap about 5 mm, for ®^Rb the dimensionless displacement will 
be = 8 X 10^ 
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After n/2 oscillations, the atoms scattered into the first order will enclose an area Tmd/uj\\ 
and the accumulated phase shift is 

= VL. (35) 

To measure the rotation rate of the Earth, with a vr phase shift in a trap with uj\\ = 27r x 
3 Hz, the trap must be moved back and forth three times. The time that it takes for the 
interferometer to measure a given phase shift does not depend on the parallel trap frequency. 
For example the time that it takes to measure a vr phase shift is 

To measure Earth's rotation, with a tt phase shift, the interferometer must have a cycle 
time of about one second. To measure a given rotation frequency, the bandwidth of the 
interferometer can only be increased by increasing the distance that the atoms are displaced 
d. For the remainder of this paper only the interference signal will be discussed for the case 
where the phase shift is zero and when the trap is not moved in the y-direction. 

Figure [2] shows the interferometer signal S", which is proportional to amplitude of the 
back scattered light, as a function of time. The sold line in Fig. [2] is Eq. fl33|) for times close 
to the first oscillation period t = t — 2n7i/uj^\ where n = 1, and with the parameters a = 0.2, 
uj^ = 10-^ = 10-2, ip = 0^ N = 7 X 10^ H = 1 and T = 30. The dots are the result of 
the DSMC, with each super particle representing 10 atoms and the signal averaged over 64 
separate runs of the DSMC code. The error bars in all DSMC calculations are smaller than 
the size of the dots shown in the figures. This figure demonstrates good agreement between 
the analytic result and our DSMC code. 

The shape of this signal illustrates the time and position varying amplitude of the density 
modulation echo relative to the probe laser. For times slightly less than one trap period 
r < 0, the nodes of the density modulation are located at the anti-nodes of the standing 
wave laser field. At precisely one trap period, the density modulation vanishes and the 
cloud returns to its initial density distribution. For times slightly larger than one trap 
period r > 0, the nodes of the density modulation are located at the nodes of the standing 
wave laser field. 

For weak pulses S < 1, the interference signal Eq. (l33l) can approximately written as 

^= -Are-^"'/'-2"^"/"ii, (37) 



13 



where A = N JkJk+i{k + l). The two peaks in the signal occurs at the times r = ±T~^/^. 
For ^^Rb at 40 fiK the time between the maximum and minimum signal is about 1 fis. The 
magnitude at the peaks in the signal is 

_ iV / oolaNn 1\ 

bpeak - ^7p^ exp I - 21/2^^2^ " 2 j ' ^ ' 

where Eq. ( IT6|) was used. For the remainder of this section, Eq. (138|) will be analyzed for 
several illustrative cases. Using this analysis, limits on the performance of the interferometer 
will be discussed. 

Figure [3] shows the maximum value of the interference signal as a function of perpendicular 
trapping frequency cj^, where all the other parameters are the same as in Fig. [2J The solid 
curve is found using Eq. (138|) and the dots were extracted from the DSMC calculation. 
This figure demonstrates excellent agreement between the analytic and DSMC models of 
the interferometer. The maximum interference signal is observed for small values of the 
perpendicular trapping frequencies. This is because the density of the atoms decreases as 
the atoms are confined less tightly in the perpendicular direction. 

From Fig. [3l it is clear that the optimal value of perpendicular trapping frequency is 
the smallest value such that the movement of the trap remains adiabatic. The minimum 
transverse trap frequency can be estimated by using Eq. (12T|) and To remain adiabatic, 
the ratio between the transverse and perpendicular trapping frequencies must be 

— > 39 

where d is the maximum displacement of the trap in the perpendicular direction and T is 
the temperature of the gas. 

Figure H] shows the maximum value of the interference signal as a function of the tem- 
perature T of the trapped gas. The remaining parameters are the same as Fig. [2l The solid 
curve is Eq. (l38l) and the dots were extracted from the DSMC calculation. There is still 
good agreement between the analytic and DSMC models. 

Holding all other parameters constant, the interference signal becomes smaller as the 
temperature is reduced. This is because, in a harmonic potential, collision rate is inversely 
proportional to temperature. For the parameters used in Fig. HI the signal increases with 
temperature until T = 60. For temperatures larger than T > 60, the duration of the echo 
becomes shorter and the amplitude of the density modulation is reduced. Using Eq. (I38l) it 
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can be shown that the largest amphtude of back-scattered hght occurs when the temperature 
is 

T = ^ . (40) 

TT 

Equation (HOl) shows that the optimal temperature increases linearly with atom number. As 
the atom number increases, the signal to noise ratio of the detected signal decreases. The 
time between the maximum and minimum amplitude decreases. The speed of the detection 
scheme places an upper limit on the temperature and therefore the lower limit on the signal 
to noise ratio. Analysis of the details of the detection scheme are beyond the scope of this 
paper and will be left to future work. 

The initial temperature of the atomic gas depends on the details of the laser cooling and 
loading of the gas into the trap. Although it is possible to experimentally vary the final 
temperature, it is easier to vary the number of trapped atoms. This can be done by changing 
the load time of the magneto optical trap. Because of this, we believe that it is most useful to 
treat temperature T, and trap frequencies u\i and uj± as constants and optimize the number 
of trapped atoms A^. 

Figure [5] shows the maximum value of the interference as a function of number of trapped 
atoms A^. The solid curve was found using Eq. (138!) and the dots were extracted from our 
DSMC code. The remaining parameters are the same as in Fig. [2l When the number of 
atoms in the gas is A^ < 7 x 10^, the signal increases with increasing number. This is because 
as the number of atoms increases so does the amount of back scattered light. When the 
total number of atoms in the gas is A^ > 7 x 10^ the interference signal decreases because 
the higher density increases the collision rate between the atoms in the gas. Using Eq. ( 138|1 
it can be shown that, holding all other parameters constant, the maximum value of the 
interference signal occurs the number of atoms is 

N = (41) 

For a trap with frequencies uj\\ = 27i x 3Hz and uj± = 2-71 x 300Hz, that traps ^'^Rb atoms at 
40 fiK the optimal number of atoms for one trap period is about 7 x lO^atoms. To measure 
Earth's rotation with at tt phase shift, by displacing the trap by 5 mm, the atoms must 
oscillate three time in this trap and the optimal number of atoms is 2.2 x 10^. 



15 



IV. CONCLUSIONS 



In this paper we presented a simple analytic model of the dynamics of a trapped atom 
interferometer that uses a single Kapitza-Dirac pulse to modulate the atoms and the classical 
turning points of the trap to reflect them. The interferometers signal is read out by reflecting 
a single probe pulse off of the atoms and interfering the back-reflected light with a reference 
beam. We presented a description of the collisions between the atoms and showed that our 
simple model give quantitatively accurate results when compared to a DSMC model of the 
interferometer. Finally, we used our model to find the optimal temperature or number to 
maximize the performance of the interferometer. 

Although the analytic model presented in this paper specialized to the analysis of a single 

ffitza-Dirac pulse, the results of Sec. [TTll easily generalize to multi-pulse interferometers 
Q. TO appl, o.. .node, to ,„t.fe.o,nete. that use gases above the BEC transition 
temperature but below the recoil temperature the momentum and spacial dependence on 
the coUision frequency cannot be ignored and Eq. (fT2l) must be used instead of Eq. (fT6l) . 
We believe that inclusion of the more complicated collision frequency will not dramatically 
change the results of this paper when describing a gas below the recoil temperature. 
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FIG. 1: (color online) A schematic of the interferometer cycle, (a) The cold gas is loaded into a 
cigar shaped trap, created with the upper most horizontal wire (b) The atomic cloud is illuminated 
with a standing wave laser field, (c) After the laser pulse, a density modulation appears across the 
cloud, (d) The density modulation disappears and the trap is moved downwards by cycling the 
current in the wires, (e) The trap is above the bottom wire at a half period of the parallel trap and 
the density modulation reappears across the cloud, (f) The trap is moved upward and the density 
modulation disappears, (g) The trap reaches the top wire at one trap period of the parallel trap 
and the density modulation reappears for a second time, (h) A probe pulse is reflected off of the 
cloud. 
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FIG. 2: (color online) The interference signal as a function of time near one trap period. The solid 
curve is Eq. ([33|) and the dots are the results of the DSMC code. 
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FIG. 3: (color online) The maximum value of the interference signal as a function of perpendicular 
trapping frequency. The solid curve is found using Eq. ()38l) and the dots were extracted from a 
DSMC calculation. 
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FIG. 4: (color online) The maximum value of the interference signal as a function of temperature. 
The solid curve is Eq. ()38p and the dots were extracted from a DSMC calculation. 
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FIG. 5: (color online) The maximum value of the interference signal as a function of number 
of atoms in the trap. The solid curve is Eq. (1381) and the dots were extracted from a DSMC 
calculation. 
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